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Abstract. We examine the possibility that the strong heating produced at 
temperature-mlnlnui:. ' v'els during solar flares Is due to resistive dlsslpa- 
tlon of Alfven waves generated by the primary energy release process In the 
corona. It Is shown how, for suitable parameters, these waves can carry 
their energy essentially undamped Into the temperature-minimum layers and 
can then produce a degree of heating consistent with observations. 
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1. Introduction 


There now exists substsntial observstlonsl evidence for strong heating 

of the temperature>mlnlmum (T . ) region In solar flares (Machado et al.. 

mxn ” 

1978; Cook» 1979; Cook and Brueckner, 1979; Canfield et al, . 1980). This 

observed temperature rise amounts to some 100-20C K, corresponding to a 

—3 —1 

Specific energy input of some 10>20 erg cm s at these levels (Machado 
et al .. 1978). It appears to be difficult to reconcile this amount of heat- 
ing tflth a primary flare energy release in the corona and subsequent transport 

» 

of this energy through the atmosphere by either particles or radiation 
(Machado et al . . 1978; BoasUe and Machado, 1979). Emslle and Machado (1979), 
among others (e.g., Sturrock, 1980), have therefore proposed that the T^^ 
region is heated by joule dissipation of currents. They do not consider 
the generation of these currents, but they do point out that such currents 
must have a very fine spatial structure 10^ cm) since otherwise the 
magnetic field associated with them would be unacceptably large. 

According to Machado et al . (1978), the total energy dissipated at T^^ 
levels during flares is small compared to the total flare energy budget. Thus, 
in explaining heating, the total energy release requirement is not 

difficult to meet. The problem Instead is to find a mechanism for transport- 
ing this energy efficiently from the primary release site to the T ^^^ level 
of the atmosphere (or to explain how part of the energy can be released at the 

Emslle and Machado (1979) note that, because of the relatively low tempera- 
tures and low fractional ionization level of the plasma at levels, the 

plasma resistivity r\ there is much larger than at any other point in the 
solar atmosphere (see Figure 1). Hence the energy carried by electric currents 
generated in the corona may reach the T^^^ layers of the atmosphere, essen- 
tially tmattenuated, and then be efflclcmtly damped due to the large resistivity. 
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The transverse scale of such coronal current patterns generated during 
flares may well be small enough (see Section 3) to satisfy the constraints 
Imposed by Smslle and Machado (1979). 

In this paper we therefore investigate the possibility that electric 
currents associated with AlfWIn^ave perturbations, generated by the primary 
flare energy release in the corona, are responsible for heating tha teiiq>era> 
tureonlnlmum region during flares. In Section 2 we discuss the propagation 
and d«*«np<ns of these waves. In Section 3 we investigate the heating rates 
effected by this mechanism for a range of model parameters and so place 
constraints on these parameters if the mechanism is to be effective at 
heating. In Section 4 we discuss the results obtained. 


Propagation and Damping of Alfven Waves in the Solar Atmosphere 
We wish to analyze the transmission of enetrgy along magnetic field lines 
from coronal heights down to the temperature -minimum region. In the MHD 
approximation, there are three modes of propagation: the slow and the fast 

mangetoacoustic modes, and the Alfven mode. In a con^lex situation involving 
non-uniform plasma and magnetic field, the group velocity of the magneto- 
acoustic modes must be expected in the course of propagation to depart from 
the direction of the magnetic field, even if the group velocity is initially 
directed along the field (Weinberg, 1962). On the other hand, the group 
velocity of Alfven waves is always parallel to the magnetic field. Energy in 
the Alfven mode will therefore remain in that mode as long as the wavelength 
is small compared with the length scale for variations of the field and plasma. 

The scale height for variation of the plas. a properties varies over a 
wide range in the sun's atmosphere. The wavelength of a propagating Alfven 
wave also varies with height since the Alfven speed varies. We find that, for 
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waves of incereait Che ws'/elengch Is much less than che scale height for 
variation of density (for Instance) throughout the corona and in the upper part of 
the transition region (see Section 3). We begin by considering this region 
and postpone until later in this section discussion of the lower transition 
region In which the density varies abruptly* 

Since we are concerned only with the propagation of energy along magnetic 
field lines* we shall ignore che magnetoacouatlc modes and consider only the 
Alfven mode. We begin by considering a simple modal in which the magnetic 
field Is uniform and vertical, ^ shall allow for non-zero reslsltlvity 
and for variation of density with height. We shall subsequently consider the 
additional effect of slow variation of the magnetic field strength with height. 

The equation of motion of the plasma Is 

p dt * ’ 

where we use modified Gaussian units (Jackson* 1962) for electromagnetic 

quantities* and we have ignored the gravitational force since we shall be 

concerned with wave frequencies much higher than the Brunt-Valsala frequency 

(Unno et al . * 1979) . The unperturbed state is characterized by the density 

p (z) , the pressure p (z) and magnetic field strength B (z) > (0,0, B ) . 
o o 'VO o 

To first order, the perturbation In the equilibrium state of the plasma and 
field satisfies 

dfiy 

“ dT • - + 

On using the Maxwell equations 
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and 


« 

and the Ohm 'a law relatlonahlp 

+ «- 5 ) 

where the resistivity n ia taken to be a scalar, we obtain 

^ * TS ^ - IS ' • «•« 

On combining equations (2.2), (2.4) and the appropriate perturbed form of 
(2.6), we obtain the wave equation 

5 ■ ^ ^ " «o]" «oj 

4? " - -JiF 

Since the Alfvin node is transverse, we may without loss of generality consider 
that ■ (0, 6B^, 0) . Further, since the Alfv6n mode involves divergence- 
free motion, the pressure gradient plays no role. (Our assumption that the 
wavelength is small compared with the scale height in fact makes it possible 
to neglect this term not only for the Alfv6n mode but also for the fast and 
slow magnetosonic modes.) Hence we may write 

- b(z)e^^*‘x* “ . (2.8) 


S 


Equation (2.7) then reduces to 


- tt>^b 


Ivl « ISSil in « icn<n 
de V dz / 4 it dz dz Air 




(2,9) 


where V is Che Alfv&t speed given by 


If (..«aice V) is 

varies as e^z^. where k 

z 

waves : 



Airp^ 

independent of z and if n ~ 0 , then b 
satisfies the dispersion relation for Alfvdn 


( 2 . 10 ) 


0)2 - V2 R2 
z 


( 2 . 11 ) 


Since, for the problem sf interest, n 0 , the Alfvdn waves will be damped. 
If V varies slowly with z , then an Altv&i wave will retain its identity 
but its amplitude will vary in the course of propagation. We may consider 
both of these effects at the same time by expressing b(z) in the form 


b(z) 


A(z) exp 



z 


k 

z 


(z')dz 



( 2 . 12 ) 


Wa find that, to lowest order in the effects of n and the variation of V , 
the expression given by (2.12) Indeed satisfies the wave equation (2.9) if 
the amplitude A(z) varies as follows: 


A(z) « A(0) 



(2.13) 


In deriving this expression, we have assumed that the wave is propagating in 
the z-direction. The "damping length" is expressible as 
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C.14) 


h 


Lx L;, 


where Lx and L„ , the “traneverae" and "parallel" damping lengtha* are 
given by the expressions 


T AttV 4irV 

"ST" 


(2.15) 


and 


L.i 



(2.16) 


and we write 1 for (k^) ^ . 

If the energy flux being carried by an Alfv^n wave Is represented by 
S , and if the actual perturbation in the magnetic field due to the wave Is 
given by the real part of expression (2.12), then 


S 



(2.17) 


Hence, from Equation (2.13), 

SCz) » S(0) exp ITPyj • <2.18) 

This equation will be used to estimate the damping of Alfvin waves in the 
neighborhood “f the temperatureHslnlmum region. 

As an AIfv4n wave of given frequency propagates down through the corona 
into the transition region, it encounters a density gradient which Is 
initially small but becomes larger and larger. At some location, the assump- 
tion that the wave length is small compared with the characteristic length 
scale of variation of density becomes invalid. Below that point, the density 
gradient becomes very steep, so that there Is in fact a sharp change in 
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density in less thsn one wavelength. In this situation, part of the wave 
energy will be transmitted and part will be reflected. 

The component of the wave vector parallel to the z-axls is given by 

k • w (4ffm n)^^^ , (2,19) 

z p 

where n Is the number density both of electrons and protons. To a fair 
approximation, the transition region may be regarded as a region In which 
the heat flux F (erg cm a ) Is constant (Athay, 1976). Then, using 
the Spltzer (1962) conductivity, the characteristic length scale for the 
temperature change , defined by 

- T (dT/dz)“^ , (2.20) 

Is given by 

Uj. - ^ (2.21) 

Where a % 10 erg cm a K ' . The gas pressure, given by 

P - 2nkT , (2.22) 

Is approximately constant In the transition region so that the density 
scale height Is equal to the temperature scale height . Hence the 
critical condition that the density gradient Is comparable with the Alfven 
wave number may be expressed as 


k^ - 1 . (2.23) 

On using Equations (2.19), (2.21) and (2.22), we find that this condition 
occurs at the "critical" layer where the density and temperature are n 

V 

and T , given by 
o 


8 


c 


(2.24) 


<nd 


T ■ 10^*^ p“V® « 

c 


(2.25) 


As far M cransalssion md raflactlon ax« concamad, cha iaportant quattcity 
la the jvmp in danaicy in going from thia critical region Into tha chromoaphara. 
which W;i write aa . 


0 .fch 

Pc “c ^ch 


(2.26) 


We aee from Equation (2.25) that thia may be expreaaed aa 


0 ■ 10^*^ P“W® b^^® « 

ch 


(2.27) 


We write the wave amplitude on the coronal aide of the tranaition region an 

i[- k_^\z-*j) - wt] 


i[k (z-z.) - wt] « X- -w 

b_ (z.t) - Ae + »B ^ 


and the wave amplitude on the chromoapheric (tranamitted) aide aa 


(2.2f) 


b^ (z.t) • Tj Ae 


(2.29) 


where the jump ia aaaumed to occur at z ■ Zj . On noting that b ia con- 
tinuous acroaa the Jump and. from Equation (2.9), that 


db db. 

«2 1 • v2 t 

- dz V dz * 


(2.30) 


m find that the reflection and tranamiasion coefficienta ere given by 
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' + I 


T .l£IL. 
® + 1 


(2.31) 


On nocing that th* enargy flux is proportl;^*^*! to b^V , w* ••• that the 
coafflclante for the reflection and tranin ilon of energy are given by 


. - - i) ' - ._4§iiL 

^ ( 9^/2 + i ) 2 » E ( qI /2 + i )2 


(2.32) 


If energy Is suddenly released at the top of a loop, then a fraction 
of that eoecgy is trananltted to the chronosphere and photosphere after a tine 
t^ , where t^ Is the time It cakes for an Mfvdn wave to propagate from the 
top of the loop to the base. After an additional time 2t^, an additional 
fraction is transmitted, etc. (Ue have assumed Chat the wave Is not 

significantly damped In the corona; this will be verified jt posteriori in Section 
3.) Hence we see that the characterlstlv decay time of the energy flow to 
Che lower atmosphere Is given by 


T 


D 



(2.33) 


For example, consider the range of values 0 10 - 100 . For this range, 

we find that ■ .27 .67 , so that ■ 1.5 - 5.0 . In a symsetrlcal 

configuration, the energy will escape equally through both feet of the loop. 

We have so far considered that Che magnetic field Is uniform and vertical. 
In a real solar flux tube, neither of these conditions will be satisfied. Wein- 
berg (1962) has discussed very generally the effect upon wave propagation of 
variations in the parameters entering the wave equation. His analysis confirms 
chat for Alfv&i waves, even In a general magnetic-field configuration, energy 
propagates along magnetic f^eld lines and that Che wave amplitude varies accord- 
ing to Equation (2.13). 
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Hovavar, Ualnbar^'a aaalyaia alao ahais light cm tha variation of tha 
wava vactor during propagation* If tia conaidar propagation along a fiald Una 
vhlch ia :urvad but liaa in a plana, wa find that tha coaponanta of tha vava 
vactor along t‘^a fiald lina (kn) and nonMl to tha plana (k^) ara un- 
af fact ad by tha curvatura, but tha rata of changa of tha coa^onant k^ 
noraal to tha fiald Una acqulraa tha following contribution 


dk 

^ - 2ic k„ , (2.34) 


whara k(s) la tha curvature the fiald Una* In order to aaaaaa tha influ- 
ence of chia effect, wa nay cemaider the cckae that a wave bagina at tha top of 
a loop (a ■ 0) with k^ * 0 * Then at tha baae of the loop, idtara a ■ L , 


ki(L> 



da 


[V(a))v(LjJ 


(2*35) 


Aa a reault of the rapid Increaae in dana<.ty aa a approachea L , the donlnant 
contribution to the integral cotuaa from the laat (tha lowaat) denalty acala 
height. If, aa ia Ukely, the angular deflection yk(a)d8 of a magnetic 
field line la aaaU in tha lowaat danaity acala height, wa uay neglect tha 
change in the wave vector due to curvature of the field linaa* 

The wave vector will alao be modified during propagatior. by the convergence 
or divergence of the field linea. If the convergence ia "one-dimanaional" in 
Che aenae Chat there ia convergence or dlt^rgence in erne of the tranaverae 
directiona (x , aay) but not in the other direction (y) , thei^ k^ acalea 
with thct magnetic field atrength but k^ ia unaffected. If, aa wtothe*’ estample, 
the field configuration haa cylindrical a'.'retetry* then both of the tranaverae 
componenta of the wave vector acale aa * We ahall, for definiteneaa, 

conaider the latter caae in the examplea to be diacuaaed in the next aection. 
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3. Modrl ParoMtara 


Table I 


Figure 1 
Figure 2 


The variation of b(z) td.ch z , aa glvan bjr ^uaclona (2.13), (2*28) 

and (2.29)« ajpeada not only on tha val' na of tha wava paraaatara <u and 

t at tha top of tha loop), but alao on tha faatctaa of tha flara loop 

itaali* Wa hava Invaatigatad four potaibla loop aodala, ahoun in Tabla I* 

Thaaa modala differ only in thair nagnatic field atrangtha and gaonatriaa, 

with the coronal field atrengtna in tha range 100^300 G and tha field atrangtha 

near tha footpointa of tha loop in tha range 100-900 G. We conaidar a narrow 

flux tube with "cylindrical" rynDatry. For tha purpoaea of calculation, a 

convaniant fom for tha variation of B with z ia tha nodal conaidarad by 

o 

Antxochoa and Sturrock (1976), 

- aec2 (f arc coa , (3,1) 

o 

(although this nodal was derived for a field configuration of aonewhat differ- 
ent synontry). The area factor ("conpresaion ratio") T for each nodal con- 
sidered ia tabulated in Tabla I. Below z » L , the backgroiuid field strength 
was taken to be uniform and equal to B^(L) . The chronospherlc variations of 
remperature, density, ionization level and recistivlty with height were 
taken fron nodal *F1* ox Machado at al . (1980) and are shown in Figure 1, 
plotted against tha coluon number density M (coT^) as independent variable, 

Tha resistivity was asaunsd classical (Spltzer, 1962), with allowance nada for 
the affect of Ion-neutral collisions near the tenperature minimum (Ensile and 
Machado, 1979). The possibility of anomalous resistivity in the region of 
flare energy release will be considered below. 

Ue solved the relevant equations for each atnospherlc nodel for a range 
of both v(a> <k)/2r) and , the resulting behavior of the quantity Q (local 
energy daposltion rate per "nit voluae) is shown in Figure 2. Ihe results in 
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the figure correspond to an initial perturbation amplitude b(0) ■ IG ; for 
other values of b(0) , the energy deposition rate scales as [b(0)]^ « The 

shaded areas in the figure represent the approximate location of the strong 
^min enhancements reported by Kachado et al . (1978); the maximum 

of Q(H) must lie in this region if the heating mechanism proposed in the 
present paper is to be effective. For values much larger than 10^ cm » 
the L„ term in Equation (2.14) completely dominates [cf. Equations (2.15) and 
(2. 16)1, and so the results are similar to those for • 10^ cm . Similarly, 
for u « 1 Hz , the term dominates, so that results for such low frequen- 
cies need not be shown either. 

It is immediately obvious from the figure that the proposed mechanism is 
effective at heating layers only if the parameters v and fall 

within the ranges 


V ^ 10 Hz 0.2) 

and 

^ 105 cm . (3.3) 

Further, only for models with high photospheric magnetic field strengths [result- 
ing in higher Alfv^n velocities and so lower damping rates in the chromosphere 
- see Equations (2.15) and (2*16)] is a satisfactory situation found (e.g. 
models c, d in the case • 10^ cm , \> ■ 1 Hz) . This result is encouraging, 

since the results of Tarbell et al . (1979) indicate that such high values of 
the magnetic field strength at the photosphere are indeed appropriate. 

The required values of v and may be compared with independent esti- 
mates of their likely values. Recent analyses of the time structure of flare- 
associated hard X-ray bursts (Dennis et al .. 1981) reveal structures with time- 
scales of the order of a few tenths of a second, with no significant structure 


13 



on cime-scaleB smaller than this. Thus, based on this data, we would expect the 
highest frequency of any flare-aaaoclated disturbance to be of order 10 Hs , 
in good agreement with the values required (Figure 2). Further, on the basis 
that magnetic reconnection in an unstable arch is the source of the primary energy 
release, we should expect to be comparable to the width of a magnetic 
"island" (Spicer, 1976, 1977). From Spicer's (1976) Equations (A. 35) and 
(4.54) (see also Furth et al., 1963), we find that 



where K is the wavenumber of the perturbation Inducing the magnetic reconnec- 
tion. Inserting numerical values and noting that K % L ^ , we obtain, using 

the coronal model given in Table I, 

^1 % 10^*^ cm (3.5) 

(see also Van Uoven, 1976). In deriving this value we have used the classical 
plasma resistivity (Spitzer, 1962). The presence of turbulence in the reconnec- 
tion region will increase n and so . For instance, is increased to 

% 1Q5 cm if the turbulent resistivity exceeds the Spitzer resistivity by a 
factor of about 2x10'* . Hence the required value 10^ cm is not incom- 

patible with the "island thickness" likely to arise during reconnection. 

As mentioned in Section 1, the observations and empirical model-fitting of 
Machado et al . (1978) indicate flare-associated temperature rises of some 100- 

200 K at T . levels, corresponding Isf. their Equations (18) and (19)1 to 
min 

a specific energy deposition rate of some 10-20 erg cm s , depending on 
whether the observed structure is "average" or "bright". From the results of 
Figure 2 (which, we recall, are norjiolized to unit initial perturbation ampli- 
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cude b(0)) , we therefore deduce that Che required value of the magnetic- 
field fluctuation generated by reconnection should have the value 

b(0) ^ SG (3.6) 

In order for Q to agree with these observatlonally deduced values. This 
value of b(0) also corresponds to a bOCal (i.e. depth-integrated) energy 
deposition rate of some SxiO^ erg crn'^s*^^ in the layers (see dashed 

curves in Figure 2); this is also in good agreement with Machado et al . *s 
(1978) estimates. 

Such a value of b(0) results in a value of b of 15 G Just below 
the transition zone [see Equation (2.31) ]. The value of b further Increases 
to % 40 G at a colum depth N % 3x10^^ cm ^ , because the increase of p and 
Che corresponding decrease of V (Equation (2.10)], leads to an Increase of b 
[Equation (2.13)1. For larger M , b falls off due to resistive damping of 
Che waves which overcomes any Increase In b due to the effect of Increasing 
density. These values of b are much smaller than the ambient field strength 
In the models of Interest, as Is necessary to Justify the linear analysis of 
Section 2. Finally, we note that, for the values of v and emerging from 
our study, the wavelength in the corona Is much less than the temperature length 
scale in that region. Also, both Lh and [Equations (2.15) and (2.16)] are, 
for classical resistivity, much larger chan L , so that the wave damping in the 
corona (z < Zj) is negligible. In fact, even if the resistivity in the region 
of magnetic reconnection exceeds the Spitzer (1962) value by the factor of ^ IQ** 
required to reconcile with the thickness of Che reconnection region (see 

above) , it is still easily demonstrated that both L,, and ^ are large compared 
Co the size of such a region. We thus see that the modeling assumptions of 
X « L_ and L_. » L made in Section 2 are Justified. 
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4. Discussion and Conclusions 

As pointed out above (Section 1)» the total energy liberated at 
level'i during solar flares Is a relatively small factor In the total energy 
budget of a flare, and therefore the constraint that the observed strong heat- 
ing at levels places on a candidate heating mechanism does not hinge 

upon the total energy supply but simply upon the existence of a mechanism to 
transport energy throughout the solar atmosphere from the corona until it is 
efficiently damped at levels. The mechanism proposed in the present 

paper permits virtually undamped propagation of energy through the corona and 
upper chromosphere, coupled vith an efficient damping mechanism at 
levels. 

We conclude from the results of the previous section that Alfv^n waves 

of moderate amplitude can, through resistive dissipation of their associated 

alternating currents, account for the observed temperature enhancements at 

T . levels during solar flares. However, it is Interesting to note that 
nla 

the propagation of energy from the primary release site In the corona to the 

T . level is relatively slow, with an associated timescale 
min 

T./%y . (4.1) 


For the coronal and chromospheric models adopted In this paper (Table 1), we 
can readily confute this integral numerically; we find that (with T measured 
In seconds) 


T 


®o,ph 


(4.2) 


where (gauss) is the strength of the photospheric magnetic field Ithe exact 

field geometry in the corona is unimportant, since the disturbance travel time 
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in the corona represents only a small contribution to T; see Equation (2.33)]. 

This is much larger than the tlme>scale quoted by Emslle et al . (1978) for an 

MHO disturbance to propagate In the upper chromosphere » because of the much 

larger densities (and correspondingly low Alfv^n velocities) appropriate to 

the T . region. Nevertheless « for the plausible range of values . 
fluxn 

^ 300 G t the timescale is much smaller than the 10 minute delay time between 
the impulsive phase of the flare and the enhancements observed by 

Machado et al . (1978). This therefore implies that, for the proposed heating 
model to be valid, magnetic disturbances must be generated in the corona (or 
elsewhere) long after the impulsive phase of the flare. It is possible that 
ttu'se disturbances can reveal themselves in other observational signatures, 
such as radio emission; however, it is unlikely that their associated currents 
will produce observatlonally detectable effects in locations other than the 
T^n ves^on, because of the very low resistivity in all other parts of the 
atmosphere . 

Finally, we note that Equation (4.2) represents a strong observational 

test of the present model. If Equation (4.2) is valid, no enliancement at the 

T , level should be detected until a time measurably later than the impulsive 
min 

phase of the flare — as indicated, for example, by the hard X~ray burst. It 
appears that data do not yet exist against which to check this consequence 
of the model. 
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Figure Cape ions 


Figure 1. Chromospheric model adopted (from Machado et al. . 1980, model *F1'). 

Shown, as a function of colunn depth M , are x the Ionization 
fraction, n - the number density, T - the temperature, h - the 
height above the photosphere (‘^3 oqq * !)• ^ ** plasma 

resistivity, assumed classical (see BmsUe and Machado, 1979). 

Figure 2. Energy deposition rate Q (solid) and depth- Integrated deposition 

rate fota (dashed) versus N . Each diagram refers to a different 

set of (v, 2. ) values as shown around the border, w'nlle within 
o 

each diagram the letters a through d refer to the coronal magnetic 
field geometry of the model (see Table 1 ). The shaded areas 
correspond to the location of the temperature minimum (cf. Figure 
1); Q should peak In this depth Interval for the model to be 
effective at heating levels of the atmosphere. 
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